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ABSTRACT 


The Milagro gamma-ray observatory employs a water Cherenkov detector to 
observe extensive air showers produced by high energy particles interacting in 
the Earth’s atmosphere. Milagro has a wide field of view and high duty cycle, 
monitoring the northern sky almost continuously in the 100 GeV to 100 TeV 
energy range. Milagro is, thus, uniquely capable of searching for very high- 
energy emission from gamma-ray bursts (GRBs) during the prompt emission 
phase. Detection of >100 GeV counterparts would place powerful constraints 
on GRB mechanisms. Twenty-five satellite-triggered GRBs occurred within the 
field of view of Milagro between January 2000 and December 2001. We have 
searched for counterparts to these GRBs and found no significant emission from 
any of the burst positions. Due to the absorption of high-energy gamma rays by 
the extragalactic background light, detections are only expected to be possible 
for redshifts less than ~0.5. Three of the GRBs studied have measured redshifts. 
GRB 010921 has a redshift low enough (0.45) to allow an upper limit on the 
fiuence to place an observational constraint on potential GRB models. 

Subject headings: gamma rays: bursts — gamma rays: observations 


1. Introduction 

The search for prompt very high energy (VHE) emission (>100 GeV) from gamma-ray 
bursts (GRBs) is motivated by both experimental observations and theoretical predictions, 
and its detection could allow us to constrain GRB emission models. Although no observation 
has yet conclusively demonstrated VHE emission from any single burst, there have been 
several indications of emission at these high energies. 
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Milagrito, a prototype of Milagro, reported evidence for emission above 650 GeV 
from GRB 970417a, with a (post-trials) probability of 1.5x10“^ of being a background 
fluctuation (Atkins et al. 2000a, 2003a). This search included 53 other bursts from which 
no signihcant emission was detected. The Tibet air shower array reported a correlation 
between ~10 TeV air showers and a sample of 57 GRBs detected by the Burst And 
Transient Source Experiment (BATSE) (Amenomori et al. 1996), although no signihcant 
signal was detected from any of the directions of the individual BATSE bursts. Evidence at 
about the 3 sigma level with HEGRA has been published for emission above 20 TeV from 
GRB 920925c (Padilla et al. 1998). Follow-up observations above 250 GeV by the Whipple 
atmospheric Gherenkov telescope (Gonnaughton et al. 1997) did not hud any high energy 
afterglow for 9 bursts studied. Because this search involved slewing the Whipple telescope 
into position after receiving a burst alert from BATSE (for these 9 bursts the delays ranged 
from 2 minutes to 56 minutes), only delayed emission could be detected. In addition, the 
small held of view of Whipple (3°), compounded with the relatively large uncertainty in the 
BATSE location of the bursts (~ 10° in diameter) meant the process of searching the burst 
region required multiple pointings, taking a 3 hour observing cycle and still covering less 
than 50% of the actual BATSE error box, further hampering any early detection of TeV 
emission (Gonnaughton et al. 1997). 

Although most GRBs have been detected in the energy range between 20 keV and 1 
MeV, a few bursts have been observed at energies above 100 MeV by EGRET, including 
the detection of an 18 GeV photon from GRB 940217, over 90 minutes after the start of 
the burst (Hurley et al. 1994), indicating both that the spectra of some GRBs extend to at 
least GeV energies and that this emission may be delayed (Dingus 1995, 2001). 

Recently, a second spectral component was found in GRB 941017 (Gonzalez et al. 
2003). This second component extended up to at least 200 MeV, and its flux decayed more 
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slowly than the lower energy component. It is still unknown to how high in energy the 
second component extends. It is not clear whether it is similar to the high energy peak seen 
in TeV sources and attributed to inverse Compton emission or the result of a completely 
different mechanism. The second component has a very hard, power law spectrum, with a 
differential photon index of -1±0.3, which if extrapolated to 100 GeV would make the burst 
extremely bright, with a fluence greater than 10“^ erg s“^ cm“^. 

The observation of VHE emission from GRBs is hindered by collisions between 
the gamma rays and the extragalactic infrared background light (EBL), producing 
electron-positron pairs (Nikishov 1961). The degree of gamma-ray extinction from this 
effect is uncertain, because the amount of EBL is not well known. Direct measurements 
of the extragalactic background light have proven difficult, because of the foreground 
contribution from the Galaxy. There are several models of the extinction (Primack et ah 
1999; Stecker & de Jager 1998; De Jager & Stecker 2002), which are similar in their general 
features because of the constraints from the available data. Recent progress in the held has 
come about due to the accurate determination of the cosmological parameters, as well as a 
greatly increased knowledge of the luminosity function of galaxies. The most recent model 
now predicts a somewhat smaller absorption than was previously expected (Primack et al. 
2004), with an optical depth predicted to be roughly unity to 500 GeV (10 TeV) gamma 
rays from a redshift of 0.2 (0.05). 

While VHE emission from GRBs has been elusive, it is a natural byproduct of most 
GRB production models and is often predicted to have a huence comparable to that at 
MeV energies (Dermer, Ghiang & Mitman 2000; Pilla & Loeb 1998; Zhang & Meszaros 
2001). This is a result of the fact that the MeV emission from GRBs is likely synchrotron 
radiation produced by energetic electrons within the strong magnetic held of a jet with bulk 
Lorentz factors exceeding 100. In such an environment, the inverse Gompton mechanism for 
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transferring energy from electrons to gamma rays is likely to produce a second higher energy 
component of GRB emission with fluence possibly peaked at 1 TeV or above. The relative 
strengths of the synchrotron and inverse Compton emission depend on the environments of 
the particle acceleration and the gamma ray production. 

In this paper, we use the Milagro gamma-ray observatory to search for VHE emission 
during the prompt emission of GRBs during 2000 and 2001. The GRBs have been detected 
by one or more of several satellite instruments: BATSE, Rossi X-ray Transient Explorer 
(RXTE), BeppoSax, High Energy Transient Explorer (HETE), and the Third Interplanetary 
Network (IPN), and occurred overhead within the held of view of Milagro. The combination 
of large held of view and high duty cycle make Milagro the best instrument available for 
conducting this type of search. 


2. The Milagro Observatory 

Milagro is a TeV gamma-ray detector which uses the water Cherenkov technique to 
detect extensive air showers produced by very high-energy gamma rays as they traverse the 
Earth’s atmosphere (Atkins et al. 2000b). Milagro is located in the Jemez Mountains of 
northern New Mexico (35.9° N, 106.7° W) at an altitude of 2630 m above sea level, and has 
a held of view of ~2 sr and a duty cycle of over 90%, making it an ideal all-sky monitor of 
transient phenomena at very high energies, such as GRBs. The ehective area and energy 
threshold of Milagro are a function of zenith angle, due to the increased atmospheric 
overburden at larger zenith angles, which tends to attenuate the particles in the air shower 
before they reach the ground. During the period covered by these observations, the Milagro 
trigger required that approximately 55 tubes be hit within a 200 ns time window. For more 
details on Milagro see Atkins et al. (2003b). Figure 1 shows the effective area to gamma 
rays of Milagro during the years 2000 and 2001. The plot shows the effective area for four 



different ranges of zenith angles: 0°-15°,15°-30°, 30°-40°, and 40°-45°, illustrating the drop 
suffered in effective area at energies below a few TeV for large zenith angles. As can be 
seen, the sensitivity of Milagro varies slowly with angle out to around 30 degrees and then 
drops off. The effective area of Milagro ranges from 10 at 100 GeV to ~10® at 
10 TeV. Systematic effects result in an uncertainty in the effective area estimated to be no 
more than 15%. The angular resolution is approximately 0.7 degrees. 

The energy response of Milagro is rather broad with no clear point to dehne as an 
instrument threshold. To obtain a rough guide of the range of energies to which Milagro is 
sensitive, we consider a power law spectrum with a differential photon index, a, of -2.4. The 
energy (E 5 ) above which 95% of the triggered events from such a spectrum are obtained is 
approximately 250 GeV, the energy (E 95 ) below which 95% of the triggered events come is 
25 TeV, and the median energy is 2.5 TeV. This illustrates both the breadth of the energy 
response of Milagro, showing that the Milagro Detector is sensitive to a low energy (< 500 
GeV) signal. 


3. The GRB sample 

Table 1 shows a summary of the sample of satellite-triggered GRBs within the held of 
view (up to zenith angles of 45°) of Milagro during 2000 and 2001. There were three bursts 
during this two year period which were within the Milagro held of view but occurred when 
Milagro was not operating. The GRB sample presented here represents approximately 
a third of all well-localized bursts during that period (Greiner 2004). Of the sample, 
three have known redshift, with GRB 010921 (z=0.45) being the closest. The detector 
conhguration was upgraded starting in 2002 and later bursts will be the topic of a future 


paper. 
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The first column in Table 1 is the GRB name, which, following the usual convention, 
represents the UTC date (YYMMDD) on which the burst took place. The second column 
gives the instrument(s) that detected the burst. In the case of BATSE, we list the trigger 
ID. The third column gives the time at which the burst triggered the particular instrument, 
in Truncated Julian Date (TJD)^^ followed by the UTC second of the day. Column four 
gives the coordinates (right ascension and declination, in degrees) of the burst. The fifth 
column gives the uncertainty in the measured location of the burst in degrees. A dashed line 
implies that the error was small compared to the Milagro point-spread function (PSF). For 
BATSE bursts, the error radius is the 90% systematic plus statistical uncertainty (added in 
quadrature). The duration listed in the sixth column is the T90 interval for that burst as 
reported by the respective instrument teams. When times were reported for several energy 
bands, we picked the duration in the highest energy band. The seventh column gives the 
redshift of the CRB, when known. Column eight describes what (if any) afterglows were 
detected in different energy bands: (O)ptical, (R)adio, or (X)-ray. If an afterglow was 
observed, it is marked with a check mark, if it was not detected it is marked with a cross, 
and if no observation was attempted it is marked with a dot. Column nine lists the zenith 
angle of the burst at Milagro, in degrees. We include only bursts for which the zenith angle 
was less than 45°. As is implied by Figure 1, the effective area of Milagro at zenith angles 
greater than 45° becomes negligible in the energy range where we expect CRB emission to 
be detectable (e.g. < 1 TeV). The remaining columns of the table list the Milagro results 
and are described later. First, we describe in a little more detail the most interesting CRB 
from the sample. 


^^TJD = Truncated Julian Date = JD - 2,440,000.5 = MJD - 40,000 



3.1. GRB 010921 


GRB 010921 was detected by the HETE satellite (Ricker et al. 2002), and together 
with data from the Interplanetary Network (IPN) it was localized to a 310-square-arcmin 
3-sigma error box. It had a duration (Tgo) of 24.6 seconds in the 7-400 keV energy range 
and a fluence in the 30-400 keV range of 1.02x10“^ erg cm“^ (Barraud et al. 2003). A 
variable optical source exhibiting a power law decay measured by the Large Format Camera 
on the Palomar 200-inch telescope was identihed as the afterglow of GRB 010921 (Price 
et al. 2001). This was superimposed on a bright galaxy, assumed to be the host galaxy, 
determined to be at z=0.45 (Djorgovski et al. 2001). This was the hrst afterglow detected 
from a HETE-localized burst. Observations with the Hubble Space Telescope failed to 
detect a coincident SN to a limit 1.33 magnitude fainter than SN 1998bw at the 99.7% 
conhdence level, making this one of the most sensitive searches for an underlying SN (Price 
et al. 2003). A radio afterglow was detected by the VLA (Price et al. 2001). This is the first 
burst in the held of view of Milagro known to be close enough to be potentially detectable 
above 100 GeV. 

The photon spectrum of GRB 010921 in the (7-200) keV range can be ht with a 
cut-off power law, dehned by dN/dE = AE“ exp(—E/E q) with a equal to —1.49lg;g, and 
Eq equal to 206^48 (Barraud et al. 2003). Sakamoto et al. (2004) obtained slightly 
different parameters {a = —l.SSlg;^®, Eq = 197^(34 keV) by doing a joint spectral ht 
using WXM in addition to FREGATE data. The data can also be ht by the Band 
function (Band et al. 1993), dehned by dN/dE = AE"exp(—E/Eg) for E < {a — P)Eq, and 
dN/dE = BE^ for E > {a - (3) Eq, where B = A[{a - (3) x x exp {(3-a). The 

best ht parameters obtained from this ht were a low-energy power law index a equal to 
—1.52lg;Q9, a high-energy power law index [3 equal to —2.33](7;g7, and Eq equal to lOSj^gg 
keV (T. Q. Donaghy 2005, private communication). The break energy Eq is related to the 
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peak energy in the vFjj spectrum by the equation Ep = Eq ^ {2 + a). The high-energy 
power law index, (3, is not very well constrained by the HETE data, as seen by the results 
of the £t. Figure 4 shows the shape of the Band function with (3 values equal to -1.99, 
-2.33, and -2.67. As (3 becomes increasingly negative, the high-energy power law eventually 
becomes irrelevant in the HETE energy range as the fitting function reduces to the simple 
absorbed power law. 


4. Data Analysis 

A search for an excess of events above those expected from the background was 
made for each of the 25 bursts in the sample. The search was performed for two different 
durations, T90 (the time interval over which ninety per cent of the flux is observed) and 
five minutes. Although the value of T90 is derived from observations made at much lower 
energies than Milagro detects, EGRET showed that this duration is relevant at higher 
energies too. Four GRBs which were observed with EGRET were among the hve brightest 
bursts observed by BATSE and the significance of the EGRET detections in the T90 
interval ranged from 6 to 12 sigma, leading to the speculation that all GRBs might have 
high energy emission during their resepctive T90 time intervals, and EGRET simply did 
not have the sensitivity to detect the rest of the bursts (Dingus 2001). Figure 2 shows 
the distribution of significances of the 17 well-localized bursts for each TOO duration. The 
distribution plotted can be £t with a gaussian of mean -0.3 and standard deviation of 1.1. 
We also chose to search for emission from these GRBs for a duration of five minutes. This 
timescale was motivated by the recent discovery of a second higher energy component in 
GRB 941017. While the TOO for that burst was 77 s, the second component (which has a 
fluence more than three times greater than the fluence in the BATSE energy range alone) 
has a duration of approximately 211 seconds (Gonzalez et ah 2003). 
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The total number of events falling within a circular bin of radius 1.6° at the location 
of the burst was summed for the duration of the burst and is shown in column ten (Nqn) 
of Table 1. An estimate of the number of background events was made by characterizing 
the angular distribution of the background using two hours of data surrounding the 
burst (Atkins et ah 2003b). This number is given as Nqff in column eleven of Table 1. 
The signihcance of the excess (or deficit) for each burst was evaluated using equation 
[17] of Li & Ma (1983) and is given in column twelve. For bursts with an uncertainty in 
the position greater than 0.5°, the search region (given by column hve of Table 1) was 
tiled with an array of overlapping 1.6° radius bins spaced 0.1° apart in right ascension 
and declination and the numbers given in columns ten and eleven are chosen as the ones 
giving the excess with maximum significance. Since we do not take into account the total 
number of bins searched, this number does not represent the signihcance of detection of this 
burst and we therefore do not include these values with those of the well-localized bursts, 
in column 12. No signihcant excess was found from any burst in the sample. The 99% 
conhdence upper limits on the number of signal events detected, Nul, given the observed 
Non and the predicted background Nqff, following the method described by Helene (1983) 
are given in column thirteen. Finally, we convert the upper limit on the counts into an 
upper limit on the huence. Using the effective area of Milagro, Aeff, and assuming a 
differential power-law photon spectrum we integrate in the appropriate energy range and 
solve for the normalization constant. We chose a spectrum of the form dN/dE = 
photons/TeV/m^. The power law index of 2.4 was chosen as a conservative value for the 
spectrum. The average spectrum of the four previously mentioned bright bursts observed 
by EGRET had a differential photon spectrum with index 1.95±0.25 (Dingus 2001) and 
the average spectrum of high energy blazars can be ht with a simple power law with 
an index of 2.2 (Weekes, T. C. 2003). The normalization factor K can be calculated by 
solving the following equation: Nul = J Aeff{dN/dE)dE. Finally, we integrate the photon 



13 


spectrum multiplied by the energy to obtain the corresponding value for the total fluence: 
F = j E{dN/dE)dE, integrating from 0.25 to 25 TeV. 


5. Results 

Table 1 shows a summary of GRBs between January 2000 and December 2001 analyzed 
with Milagro. None of these bursts showed any significant emission in the Milagro data. 
The 99% conhdence upper limits on the fluence for both durations are listed in columns 
fourteen and fifteen of Table 1 respectively. Note that none of the results listed in the 
table take into account the effect of absorption from the EBL. For GRB 000301G and GRB 
000926, both of which are at z~2, we expect essentially all of the TeV emission (if it exists) 
to be absorbed by the EBL, making the actual upper limit on the emission of TeV gamma 
rays from these bursts considerably weaker than what is reflected in the table. For all the 
other bursts (except GRB 010921) there are no redshift measurements, so it is not possible 
to take this effect into account. The analysis of GRB 010921, which we have described in 
detail, yields an interesting upper limit, as we describe in the next section. 


5.1. GRB 010921 

GRB 010921, with a measured redshift of z=0.45, is close enough that a potential 
TeV signal is not totally absorbed by the extragalactic background. Figure 3 shows the 
effect of various absorption models at a redshift of 0.45 on an E“^'^ spectrum which might 
be expected from a GRB. For an unbroken power law, the energy range from 0.25 TeV 
to 25 TeV is the sensitive region for Milagro. For a cutoff spectrum, the small amount of 
sensitivity between 50 and 250 GeV is not negligible, so we now include it. Integrating the 
product of a normalization constant K times an E“^'^ energy spectrum multiplied by the 
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EBL absorption and effective area in the energy range 0.05-100 TeV and setting it equal to 
the 99% conhdence upper limit on the counts (20.2), allows us to solve for K (2.0x10“^ 
photons TeV“^ m“^). We then use this normalization constant to compute the unabsorbed 
fluence. The upper limits obtained are < 2.9 x 10“® erg cm“^ for the Primack et ah (2004) 
model, 5.8 x 10“^ erg cm“^ for the Primack et ah (1999) model, and < 5.8 x 10“^ erg 
cm“^ for the “fast evolution” model of Stecker & de Jager (1998) (the more absorptive of 
their two models). Figure 4 shows the spectrum of GRB 010921 obtained by HETE (in 
the energy range of 7 to 200 keV, where the data were htted) as well as the Milagro upper 
limits. The quantity plotted is E^dN/dE, and has been evaluated at the median energy of 
events that would be detected resulting from using each of the EBL absorption models, 
approximately 100 GeV for both the “fast evolution” model of Stecker & de Jager (1998) 
and the Primack et al. (1999) model, and 150 GeV for the Primack et ah (2004) model. 
The derived upper limits of E^dN/dE, using these three absorption models are 1.6xl0“^erg 
cm“^ (for the hrst two models) and 6.7xl0“®erg cm“^ (for third one). These are shown as 
two separate arrows in Figure 4. 


6. Discussion 

The conclusive detection of TeV emission from GRBs would improve our current 
understanding not only of GRBs themselves, but also of a whole range of physical 
phenomena. Besides possibly allowing us to learn something about the magnetic fields and 
electron energies in the GRB environments, as in the case of blazars, it could also serve to 
distinguish particle acceleration models, constrain the models of the infrared photon density 
of the universe, and possibly to probe some more exotic phenomena such as quantum 
gravity (Amelino-Gamelia et al. 1998). The number of GRBs visible to Milagro so far has 
been relatively small. The fraction of GRBs expected to be detected by Milagro depends 
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on many assumptions. As mentioned in the introduction, there are good theoretical as well 
as experimental reasons to expect GRBs to emit photons at GeV to TeV energies. Even 
if these photons were created at (and escaped) the source, however, one of the greatest 
limitations in their detection involves the attenuation due to the intergalactic background 
radiation. High energy 7 -rays interact with intergalactic star light, creating electron 
positron pairs (Primack et ah 1999; Stecker & de Jager 1998; De Jager & Stecker 2002; 
Primack et al. 2004). The attenuation is a function of gamma-ray energy, the density and 
spectrum of the background radiation helds and, crucially, the distance to the source. Even 
at the relatively modest redshift of 0.5, the optical depth is greater than 1 at a few hundred 
GeV. Figure 3 shows the effect of the EBL on an spectrum (shown, unabsorbed, as a 
straight line) at a redshift of 0.45, according to different models. Unfortunately, only three 
of the 25 GRBs studied here have a measured redshift, and of these, only GRB 010921, 
with a redshift of 0.45, was close enough for a signihcant fraction of photons to escape 
absorption. 

Figure 5 compares the Milagro upper limits for a number of bursts with the observations 
made by BATSE and HETE (we show those for which a fluence in the lower energy bands 
was available). The quantity plotted is E^dN/dE. We have assumed an E“^ spectrum at 
the lower energies and an E“^'^ spectrum for the Milagro energy range, evaluating it at 
2.5 TeV. The circles represent HETE bursts while squares represent BATSE bursts. It is 
important to note that the absorption from the EBL has not been taken into account. 
Because the redshifts of these bursts are unknown, the Milagro upper limits cannot place 
tight constraints on the burst emission. For many bursts, the Milagro fluence limits are 
below those of HETE or BATSE, so that if the burst were known to be close, the limit 
would be constraining. GRB 010921 was the hrst HETE burst with a measured afterglow. 
The VHE limit from Milagro is the hrst which constrains the TeV hux rather than the 
distance to the source. 
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With the launch of Swift, we expect ~50 well-localized bursts in the Milagro held 
of view during the hrst two years, most of them with measured redshifts. There are 
currently 39 GRBs with known redshift (Greiner 2004). Figure 6 shows a histogram of 
these redshifts (left scale), with the cumulative distribution also shown on the plot (right 
scale). Approximately 20% of the measured redshifts are 0.5 or less. A sample of known 
low redshift GRBs, combined with the upgrades made to Milagro, increasing the effective 
area of Milagro and lowering its energy threshold, should allow us in the near future to 
improve greatly on the current results and shed some light on the nature of VHE emission 
from GRBs. 


7. Conclusions 

A search for very high energy emission from GRBs was performed with the Milagro 
observatory in the range of 50 GeV to 100 TeV. A total of 25 satellite-triggered GRBs 
were within the held of view of Milagro in the two year period between January 2000 and 
December 2001, including GRB 010921, at a known redshift of 0.45. No signihcant emission 
was detected from any of these bursts. 99% conhdence upper limits on the huence are 
presented. 

Many people helped bring Milagro to fruition. In particular, we acknowledge the efforts 
of Scott DeLay, Neil Thompson and Michael Schneider. We have used GGN Notices to 
select raw data for archiving and use in this search, and we are gratefnl for the hard work 
of the GGN team, especially Scott Barthelmy. We thank Tim Donaghy and the rest of the 
HETE team for helpfnl comments as well as access to HETE data for GRBO10921. This 
work has been snpported by the National Science Eoundation (nnder grants PHY-0075326, 
-0096256, -0097315, -0206656, -0245143, -0245234, -0302000, and ATM-0002744) the US 



Department of Energy (Office of High-Energy Physics and Office of Nuclear Physics), Los 
Alamos National Laboratory, the University of California, and the Institute of Geophysics 
and Planetary Physics. 
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Fig. 1.— Effective area of Milagro to gamma rays as a function of energy for various zenith 
angles. The different lines reflect the effective area for different ranges of zenith angles (in 
decreasing order of thickness, 0°-15°,15°-30°, 30°-40°, and 40°-45°) illustrating the drop 
suffered in effective area at energies below 1 TeV for large zenith angles. 
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Fig. 3.— Effect of the EBL on an E“^'^ spectrum for z=0.45. The solid lines represent (in 
order of decreasing thickness) the unabsorbed spectrum, the model of Primack et ah (2004), 
and the Baseline and Fast Evolution models of Stecker & de Jager (1998). The dashed lines 
represent the previous four curves divided by the unabsorbed spectrum; i.e. the attenuation 
factor. 
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Fig. 4.— The GRB 010921 spectrum as measured by HETE (shown in the range 2-400 
keV). The function is the result of htting the data to the Band function. The high-energy 
power-law index f3 is unconstrained. We plot it for three values: -2.33, -1.99, and -2.67. 
The arrows are the upper limits from Milagro for various EBL absorption models, the lower 
energy one corresponding to the Stecker & de Jager (1998) and Primack et al. (1999) models 
(the results are very similar) and the slightly higher energy one corresponding to the more 
recent Primack et al. (2004) model. The energies at which these are evaluated are the 
respective median energies of the detected events and the quantity plotted is E'^dN/dE. 
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Fig. 5.— Comparison of BATSE (squares) and HETE (circles) fluences to Milagro upper 
limits. The bursts represented in this plot are: GRB 000113, 000212, 000226, 000302, 000317, 
000331, 000508, 010613, 010921, 011130, and 011212. The dotted line represents the line 


y=x. 
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Measured redshift 


Fig. 6.— Distribution of the 39 currently measured redshifts, including bursts GRB 970228 
through GRB 041006. The scale on the right refers to the cumulative distribution also shown 
in this plot. 
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Table 1. List of GRBs in the field of view of Milagro in 2000 and 2001 


GRB 

ID 2 

Time^ 

Pos> 

Err. circ.® 

T90/Dur.(s) 


ORX® 

s'* 

Non"° 


^12 

Nul^" 

UL (TgO/Dur.)"^^ 

UL (5 min)^^ 

000113 

7948 

11556:34202.4 

163.27,+ 19.89 

5.0 

370 



20.9 

471 

423.2 


95.8 

5.5e-6 

4.8e-6 

000212 

7987 

11586:81065.2 

16.09,+35.62 

4.2 

8 



2.21 

32 

17.7 


30.1 

l.le-6 

5.0e-6 

000226 

7998 

11600:13909.4 

143.68,+29.82 

4.5 

10 



31.5 

25 

11.4 


27.9 

3.4e-6 

l.le-5 

000301C 

8005R/I 

11604:35497 

245.09,+29.42 

- 

14 

2.03 


37.6 

11 

10.3 

0.2 

12.0 

2.6e-6 

6.6e-6 

000302 

8008 

11605:10225.1 

58.20,+54.28 

4.15 

120 



31.9 

143 

113.1 


54.7 

6.8e-6 

1.2e-5 

000317 

8039 

11620:77953.8 

27.22,+32.66 

5.5 

550 



6.39 

1740 

1625.9 


208.0 

7.9e-6 

5.4e-6 

000330 

8057 

11633:75449.4 

358.31,+39.26 

5.4 

0.2* 



30.0 

3 

0.3 


9.8 

l.Oe-6 

1.2e-5 

000331 

8061 

11634:85421.8 

32.00,+59.77 

5.9 

55 



38.3 

94 

66.2 


53.0 

1.2e-5 

2.1e-5 

000408 

1/8069 

11642:9348.2 

137.32,+66.58 

0.5 

2.5 


X.. 

31.1 

4 

3.4 

0.3 

8.6 

l.Oe-6 

6.7e-6 

000508 

8099 

11672:77419.3 

89.87,+2.41 

4.3 

30 



34.1 

25 

15.3 


24.0 

3.7e-6 

9.4e-6 

000615 

S 

11710:22704 

233.14,73.79 

- 

10 


XXX 

39.0 

3 

7.2 

-1.8 

6.3 

1.6e-6 

8.0e-6 

000630 

I 

11725:1853 

221.81,+41.22 

- 

20 


Vx. 

33.2 

23 

22.6 

0.1 

15.5 

2.2e-6 

6.7e-6 

000727 

I 

11752:70956 

176.00,+ 17.41 

- 

10 


X X . 

40.8 

5 

6.2 

-0.5 

8.2 

2.6e-6 

1.2e-5 

000730 

I 

11755:255 

191.29,+19.27 

- 

7 



19.2 

8 

13.9 

-1.7 

7.9 

4.2e-7 

2.7e-6 

000926 

I 

11813:85773 

256.06,+51.78 

- 

25 

2.04 


15.9 

60 

56.8 

0.4 

25.2 

1.2e-6 

3.4e-6 

001017 

I 

11834:80346 

272.18,-2.99 

0.5 

10 



42.1 

2 

4.6 

-1.3 

6.1 

2.2e-6 

l.Oe-5 

001018 

I 

11835:61114 

198.54,+ 11.81 

0.5 

31 


V- 

31.8 

31 

31.5 

-0.1 

16.9 

2.1e-6 

6.8e-6 

001019 

I 

11836:86375 

257.93,+35.34 

- 

10 


X .. 

19.5 

27 

21.1 

1.2 

20.9 

l.le-6 

l.le-6 

001105 

I 

11853:59128 

195.3,+35.49 

- 

30 


X.. 

8.5 

87 

76.2 

1.2 

35.6 

1.4e-6 

2.0e-6 

010104 

I 

11913:62489 

267.4,+18.23 

0.5 

2 



19.8 

3 

3.3 

-0.2 

7.5 

4.0e-7 

2.8e-6 

010220 

s 

11960:82267 

39.4,+61.7 

- 

150 


X X . 

27.0 

155 

168.8 

-1.1 

24.9 

2.1e-6 

3.0e-6 

010613 

I/H 

12073:27234 

255.18,+14.27 

- 

152 



24.7 

277 

280.5 

-0.2 

40.7 

2.9e-6 

5.1e-6 

010921 

I/H 

12173:18950.6 

344.0,+40.93 

- 

24.6 

0.45 

w 

10.4 

61 

65.4 

-0.6 

20.2 

8.0e-7 

2.2e-6 

011130 

H 

12243:22775.7 

46.4,+3.8 

- 

83.2 


XXX 

33.7 

93 

100.9 

-0.8 

23.0 

3.4e-6 

8.3e-6 

011212 

H 

12255:14642 

75.051,+32.13 

■ 

84.4 



33.0 

132 

113.1 

1.7 

43.8 

6.7e-6 

1.2e-5 


Note. — 2) BATSE ID, when known, otherwise: R — RXTE, 1 — IPN, H — HETE, S — BeppoSAX; 3) Time of burst, TJD: seconds; 4) ra,dec in degrees; 5) Error circle in degrees (a 
dash implies error is small compared to the Milagro PSF); 7) redshift (when known); 8) afterglow observed: O — Optical/IR, R — Radio, X — X-ray; 9) Zenith angle, in degrees; 10) Milagro 
number of counts on source; 11) Estimated number of background counts; 12) Statistical significance of the excess (in standard deviations); 13) 99% Upper Limit on the number of counts 
coming from the source; 14) 99% Upper Limit on the fiuence (0.25—25 TeV), in ergs cm~^ for the duration in column 6; 15) 99% Upper Limit on the fluence (0.25—25 TeV), in ergs cm~^ 
for a duration T = 5 min. *) Estimate. 














